Superconductivity-induced optical anomaly in an iron arsenide 
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One of the central tenets of conventional theories of superconductivity, including most models proposed for 
the recently discovered iron-pnictide superconductors, is the notion that only electronic excitations with energies 
comparable to the superconducting energy gap are affected by the transition. Here we report the results of a com- 
prehensive spectroscopic ellipsometry study of a high-quality crystal of superconducting Bao.6sKo.32Fe2As2 
that challenges this notion. We observe a superconductivity-induced suppression of an absorption band at an 
energy of 2.5 eV, two orders of magnitude above the superconducting gap energy 2A ~ 20 meV. Based on 
density-functional calculations, this band can be assigned to transitions from As-p to Fe-d orbitals crossing the 
Fermi surface. We identify a related effect at the spin-density-wave transition in parent compounds of the 122 
family. This suggests that As-p states deep below the Fermi level contribute to the formation of the supercon- 
ducting and spin-density-wave states in the iron arsenides. 

The standard Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity based solely on an effective attractive interaction 
between electrons mediated by phonons does not provide a satisfactory explanation of the properties of strongly-correlated high- 
temperature superconductors. Theoretical proposals going back many years suggest that electronic excitations might enhance 
this interaction and thus contribute to the formation of the superconducting condensate OHS). These proposals appeared to 
gain some ground with the observation of superconductivity-induced transfer of the optical spectral weight in the cuprate high- 
temperature superconductors which involves a high-energy scale extending to the visible range of the spectrum |4). In spite 
of numerous studies (for a comprehensive list of references see ref. (5) no modification of interband optical transitions in the 
superconducting state has been directly identified in the cuprates. Instead, the observed superconductivity-induced anomalies 
in the optical response of highly conducting Cu02 planes were found to be confined to the energy range corresponding to 
transitions within the conduction band below the plasma edge. These changes are dominated by a narrowing of the broad Drude 
peak caused by superconductivity-induced modification of the scattering rate (6]-[8]. Minute redistribution of the spectral weight 
between the conducting and high-energy Hubbard bands generated by Coulomb correlations may also play a role ll9l [T0l . 

Current research on the recently discovered iron-pnictide superconductors [11] suggests that electronic correlations are weaker 
than those in the cuprates. Unlike in cuprates, the Fermi surface has been reliably determined over the entire phase diagram and 
shows good agreement with density functional calculations. The superconducting state of the iron pnictides appears to fit well 
into a BCS framework in which phonons, which in these compounds interact only weakly with electrons fl2l, are replaced by 
spin fluctuations |[T3l . The ellipsometric data we present here are consistent with the hypothesis that electronic correlations 
result in only a modest renormalization of the electronic states. However, superconductivity-induced optical anomalies involve 
modification of an absorption band peaked at an energy of 2.5 eV, two orders of magnitude larger than the superconducting gap 
2A w 20 meV. In contrast to cuprate superconductors, this high-energy anomaly has a regular Lorentzian shape in both the 
real and imaginary parts of the dielectric function and is confined to energies well above the plasma edge hu p i ps 1.5 eV. It 
can be explained as a consequence of non-conservation of the total number of unoccupied states involved in the corresponding 
optical transitions due to the opening of the superconducting gap. This implies that unconventional interactions beyond the BCS 
framework must be considered in models of the superconducting pairing mechanism. 

The measurements were carried out on a single crystal of Bai_ a; K I Fe2As2 (BKFA) with x = 0.32 and superconducting T c = 
38.5 K. Specific-heat measurements on the same sample confirm its high purity and the absence of secondary electronic phases 
lfl4l . We performed direct ellipsometric measurements of the in-plane complex dielectric function e(oS) = e\{bS) + ie^w) = 
1 + 4nia(uj)/uj over a range of photon energies extending from the far infrared (Huj = 12 meV) to the ultraviolet (fax) — 6.5 eV) 
with subsequent Kramers-Kronig consistency analysis (see Supplementary Information). The far-infrared optical conductivity, 
is dominated by the opening of a superconducting gap of magnitude 2A m 20 meV below T c (Fig. [lj), in accordance with 
previous studies of optimally doped BKFA 1 15 1. The low-energy missing area in the optical conductivity spectrum below T c , 
5Al, = (af OK (oj) — al° K (uj))duj, is contained within 10A and amounts to uj^ = ^/86A[ = 0.9 eV, equivalent to a 
London penetration depth of A p = 2200 A. The fraction of the missing area below 12 meV not accessible to the experiment was 
accurately quantified from the requirement of Kramers-Kronig consistency of the independently measured real and imaginary 
parts of the dielectric function. 

Careful examination of the visible range uncovered superconductivity-induced suppression of an absorption band at 2.5 eV. 
Figure [TJ) shows difference spectra between 40 and 10 K of the real parts of the optical conductivity and dielectric function. 
The suppressed band has a Lorentzian lineshape and appears abruptly across the superconducting transition, consistently in 
both A<7i and Aei, as shown in Fig. [TJ: for <7i(2.5 eV). The temperature dependence of the suppression (blue open circles) 
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FIG. 1: a, Real part of the far-infrared optical conductivity of Bao.6sKo.32Fe2As2 and the missing area, b, Difference spectra of the real part 
of the optical conductivity (top panel) and dielectric function (bottom panel) between 40 and 10 K, with a small background shift (horizontal 
dashed line) detected by temperature modulation measurements (see Supplementary Information). Lorentzian fit to both spectra (black solid 
lines), c, Temperature scan at 2.5 eV. Contribution of the normal-state dynamics (solid line) was estimated to determine the magnitude of the 
SC-induced jump, d, Density of states in the normal (gray line), conventional superconducting state (black line), and an unconventional state 
with a depletion of unoccupied states (red line). Filled areas of respective colors represent total number of unoccupied states. 

coincides with that of the far-infrared optical conductivity due to the opening of the superconducting gap (red filled circles). Thus 
the onset of superconductivity not only modifies the low-energy quasiparticle response, but also affects the overall electronic 
structure including interband transitions in the visible range of the spectrum. Since the spectral-weight (SW) loss 5 An is not 
balanced in the vicinity of the absorption band (Fig.fTJ?), our data indicate a SW transfer over a wide energy range. We note that 
superconductivity-induced modification of the lattice parameters only results in a minute volume change of AV/V f=s 5 x 10 -7 , 
which is insufficient to explain the optical anomaly |16|. A Rramers-Kronig consistency analysis could not be carried out 
with sufficient accuracy to show whether or not the SW liberated from the absorption band at the superconducting transition 
contributes to the response of the superconducting condensate at zero energy. We did, however, detect a minute rise of the 
background level of ei (1 .5-3.5 eV) below T c , which according to the Rramers-Kronig relation implies that the SW is transferred 
to energies below 1.5 eV. This effect was identified from a simultaneous fit of A<7i(a>) and Asi(uj) (horizontal dashed line in 
the bottom panel of Fig. [TJ?). Further accurate temperature modulation measurements of o\ and E\ at characteristic energies 
confirmed a background increase of Aei(2.5 eV) = (8 ± 4) 10 -4 (see Supplementary Information). 

Since the spectral weight SAn liberated from the absorption band upon cooling below T c comprises only ~ 0.5% of the total 
spectral weight 5A^ of the superconducting condensate, its contribution to the low-energy charge dynamics might be considered 
negligible. However, assuming that this additional high-energy SW contributes to the itinerant-carrier response below T c , a 
simple estimate in the framework of the tight-binding nearest-neighbour approximation ifTTl [T8l shows that this would lead to 
a reduction of electronic kinetic energy of 0.60 meV/unit cell in the superconducting state (see Supplementary Information). 
This is close to the condensation energy AF(0) = 0.36 meV/unit cell obtained from specific-heat measurements on the same 
sample lfl4l . It is thus important to establish the origin of this unusual optical anomaly. 

We therefore compared our data to the results of ab-initio electronic structure calculations in the framework of the Local 
Density Approximation (LDA) (Figs. [2^ and[2}}). A dispersion analysis of the the experimental optical conductivity in the range 
0.5 — 6.5 eV yielded three major interband transitions in Ba .68Ko.32Fe 2 As 2 . Comparison to the LDA results enabled us to 
identify the initial and final states of these transitions. The lowest-energy transition is located at about 1 eV (red line) and stems 
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FIG. 2: a, Real part of the optical conductivity and major interband transitions of Bao.6sKo.32Fe2As2 determined by a dispersion analysis, 
b, Corresponding LDA calculation with a breakdown into separate orbital contributions described in the legend a. c, Band structure from 
the same LDA calculation. Color coding of the dispersion curves corresponds to the text color in legend a. Superconductivity-suppressed 
absorption bands (black arrows). 



from intraband Fe — d and interband As — p to Fe — d transitions. The major contribution to the optical response in the visible 
spectral range comes from transitions starting from Fe — d or As — p orbitals into strongly hybridized Fe — d to As — p or Fe — d 
orbitals (green line). Finally, the UV absorption comes from higher-energy transitions into Ba — d states (blue line). 

Although the high-energy electronic structure of BKFA is predicted quite well by the LDA calculations, the experimental 
quasiparticle response due to transitions within the conduction band (or, given the multiorbital structure of iron pnictides, a 
narrowly-spaced set of conduction bands) shows a significant deviation. The discrepancy can be quantified by the squared 
ratio of the band-structure plasma frequency ojp^ A — 2.7 eV (not included in Fig.jzj}) to its experimental counterpart uj^ p = 
1.5 eV, which can be obtained in practice from the residual optical response, after the interband transitions identified using 
a dispersion analysis have been subtracted. This ratio approximates the quasiparticle effective-mass renormalization factor 

m*/mband = ( w pi DA / w pi) ~ 3. Such an enhancement is consistent with de Haas-van Alphen and photoemission experiments 
on other compounds of the 122 family IT9142T1 and was recently reproduced by combined LDA+DMFT calculations for both 
1111 and 122 compounds El . These calculations do not show evidence of formation of Hubbard bands and thus indicate 
moderate electron-electron correlations. This explains the good agreement of the LDA optical conductivity above 1.5 eV with 
the experimental data. 

Now we turn to the physical origin of the superconductivity-suppressed absorption band. The same LDA calculation revealed 
a set of interband transitions centered at 2.5 eV, which originate or terminate in states exhibiting hole dispersion and crossing 
the Fermi level at the T- and M-points of the Brillouin zone (Fig. |2j;). We confidently assign these states to Fs-d yz ^ zx and 
Fe-d^j, orbitals. The other states involved in these transitions belong to As-p^. y /Fe-d z 2 hybrid orbitals about 2 — 3 eV below 
and above the Fermi level, giving rise to a bandwidth of AE 1 eV, in remarkable agreement with experiment. Suppression of 
the absorption band over its full width can be explained by redistribution of the occupation of Fe-d yZyZX and Fe-d xy states under 
the Fermi level below the superconducting transition. This mechanism is supported by LDA calculations in which the density of 
states within one superconducting gap energy above the Fermi level was eliminated (see Supplementary Information) leading to 
the observed suppression of the optical transitions shown in cyan in Fig.[2|). 

The required population redistribution is, however, at variance with the conventional theory of superconductivity. In the 
framework of the standard BCS approach, opening of an energy gap in a single-band superconductor leads to a bending of 
the quasiparticle dispersion and an expulsion of the density of states in the vicinity of the Fermi surface (gray and blue areas 
in Fig. [TJl) [23] and does not lead to population redistribution, i.e. the total number of unoccupied states below the transition 
is conserved = AT-JJg (blue area is equal to the gray area in Fig. [TJl). This can only lead to a small corrugation of an 
optical absorption band on the scale of one superconducting-gap energy superimposed on the overall broad feature without any 
modification of its spectral weight [24|. The experimentally observed suppression of an absorption band on the scale of its full 
width necessarily requires population imbalance A^g < N^g (red area unequal to the gray area in Fig. [TJl). This effect can be 
clearly identified as a consequence of superconductivity because the temperature dependence of the suppression mimics that of 
the optical conductivity in the FIR region due to the opening of the superconducting gap, as shown in Fig.[]J. 

All of the iron-pnictide superconductors are known to have multiple superconducting gaps [ 1 1 1 and theoretical work indicates 
a dominant contribution of electron pairing between different bands to the formation of the superconducting state [ 13 1. Redistri- 
bution of the occupation of the different bands below T c could explain the optical anomaly we observed. It requires a lowering 
of the material's chemical potential in the superconducting state. However, even a generalization of the standard BCS theory 
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to the multiband case [25 1 does not take into account this effect. Therefore, self-consistent treatment of a variable chemical 
potential at the superconducting transition is needed. In the presence of large Fe-As bond polarizability 11261 it can potentially 
enhance superconductivity in iron pnictides. We have further explored the validity of this scenario by repeating our ellipsometric 
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FIG. 3: a, Real and imaginary parts of the dielectric function of SrFe2As2- Lorentzian fit to both (solid lines), b, Temperature scan at 3.4 eV. 
c, Spectral weight (SW) redistribution between 200 K and 175 K. Extrapolation-dependent SW before (blue) and after (red) a Kramers-Kronig 
(KK) consistency check. Blue and red filled areas represent regions of SW gain and loss, respectively, in the magnetic versus normal state, 
inset, Real part of the dielectric function and KK transformations of the real part of the optical conductivity (solid lines, colors match filled 
circles). 

measurements on parent compounds of the 122 family of iron arsenide superconductors. Since the spin-density-wave (SDW) 
instability exhibited by these compounds is also believed to be induced by nesting of electronic states on different electronic 
bands, we expect a similar optical anomaly at the SDW transition as the one we observed in the superconductor. The magnitude 
of the anomaly is expected to be larger than the one in the superconductor, because the SDW transition occurs at a higher tem- 
perature and generates a larger energy gap. In SrFe2As2 (SFA), we indeed find a strong reduction of optical absorption upon 
cooling below Tguw = 200 K. The difference spectra of £i(w) and E2(oj) between 200 and 175 K show a double-peak structure 
with maxima at 2.4 and 3.4 eV (Fig. [3^). A temperature scan across Tsdw at the frequency of the second peak (Fig.JS]}) further 
confirms that this effect is induced by SDW formation. A direct comparison with the superconducting compound is complicated 
by a pronounced modification of the electronic structure due to the coincident magnetic and structural transitions. Nevertheless, 
certain information can be gained from the critical behavior of the in-plane spectral weig ht ASW(n) = £' &a x {u)duj. Fig- 
ure [3}; shows difference spectral weight of SrFe2As2 between 200 K and 175 K as a function of the upper integration limit fi. 
The change of the SW in the extrapolation region below 12 meV was accurately determined via a Kramers-Kronig consistency 
analysis (see Supplementary Information), as illustrated in the inset of Fig. [5};. With ASW(12 meV) = across the transi- 
tion (blue filled circle) the Kramers-Kronig transformation of the Ae2 (w) (blue line) deviates significantly from experimentally 
measured Asi(uj). Gradually increasing this SW brings them closer together until they finally coincide (red line) thus fixing 
ASW(12 meV) = 0.015 eV 2 (red filled circle). The higher-energy redistribution of the SW is broken down in Fig.|3j; into 
regions of SW gain (blue areas) and loss (red areas) in the SDW with respect to paramagnetic state. The SW lost due to the 
opening of the SDW gap [27] (the first red region) is partly transferred to the electronic excitations across the gap (the first blue 
region) and fully recovered by 1.5 eV. These processes are then followed by higher-energy redistribution in the region from 1.5 
to 4.0 eV involving the SW of the suppressed bands. It appears unlikely that such high-energy SW transfer could result from 
modification of the electronic structure due to a magnetic transition because effects of electronic reconstruction at the SDW 
transition are limited by 1.5 eV. A modification of the matrix elements at the structural transitions of sufficient strength cannot 
account for the observed suppression, because this would be accompanied by an even larger effect at higher energies, clearly 
absent in Fig. [3};. A redistribution of charge carriers between the SDW-coupled bands analogous to that in the superconducting 
compound provides a more natural explanation. The same physical reasons might explain the orbital polarization that breaks 
the degeneracy of Fe-d xz and Fe-d yz orbitals recently observed in the Ba-based parent of the same family by photoemission 
spectroscopy l28l . 

Interactions of electrons in different energy bands at the Fermi level may provide a common framework for an explanation 
of the optical anomalies in the spin-density-wave and superconducting compounds. It is important to note that these anomalies 
affect only a small fraction of the interband transitions, which involve initial states of As p-orbital character deep below the 
Fermi level. This indicates that these orbitals significantly influence electronic instabilities in the iron arsenides, possibly due to 
the high polarizability of As-Fe bonds. Our study points to optical spectral-weight transfer from high energies to below 1.5 eV 
induced by collective electronic instabilities. In the superconductor, it occurs at energies two orders of magnitude larger than the 
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superconducting gap energy, suggesting that electronic pairing mechanisms contribute to the formation of the superconducting 
condensate. 
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SUPPLEMENTARY INFORMATION 



Experimental 



The Bao.68Ko.32Fe2As2 single crystal was grown in zirconia crucibles sealed in quartz ampoules under argon atmosphere [$T]. 
Its chemical composition was determined by energy-dispersive X-ray spectrometry. The quality of the sample and absence of 
phase separation was confirmed by an estimated from the residual low-temperature electronic specific heat non-superconducting 
fraction of less than 2.4% [S0. From DC resistivity, magnetization and specific-heat measurement we obtained T c = 38.5 ± 
0.2 K. The sample surface was cleaved prior to every measurement. 

The experimental setup comprises three ellipsometers to cover the spectral range of 12 meV— 6.5 eV. For the range 12 meV— 
1 eV we used a home-built ellipsometer attached to a standard Fast-Fourier-Transform Bruker 66v/S FTIR interferometer. The 
FIR measurements were performed at the infrared beamline of the ANKA synchrotron light source at Karlsruhe Institute of 
Technology, Germany. For the MIR measurements we used the conventional glow-bar light source of a Bruker 66v/S FTIR. 
Finally, high-energy spectra 0.7 eV — 6.5 eV were measured with a Woollam VASE (Variable Angle Spectroscopic Ellipsometer) 
ellipsometer equipped with a UHV cold-finger cryostat operated at < 5 x 10~ 9 mbar chamber pressure. 

The inherent capacity of Woollam VASE ellipsometers to measure relative changes of the dielectric function on the order of 
10~ 2 was boosted to an unprecedented level of 10 -4 using temperature-modulation measurements (see Supplementary Infor- 
mation). The sample temperature was changed between 20 to 40 K with a period of 1800 s and later averaged over 24 periods to 
achive sufficient accuracy required to confirm a minute background increase of Aei(2.47 eV) = (8 ± 4) 10~ 4 . 

One of the strong advantages of spectroscopic ellipsometry over reflectometry is that independently obtained Aei(cj) 
and A(Ti(w) can be used in a Kramers-Rronig consistency check, in which independently obtained spectra Aeri(a;) and 
8p J °° §rzjjy dx must coincide. This additional constraint unique to ellipsometry allows one to determine with high accuracy 
the spectral weight in the extrapolation region beyond the experimentally accessible spectral range, in our case below 12 meV. 
This drastically reduces the extrapolation uncertainty and renders subsequent data analysis more robust (see Supplementary 
Information). 

The band-structure calculations were performed using a linear-muffin-tin orbital method in the atomic sphere approxima- 
tion [S3 1 within the LDA starting form the known crystal structure of Ba(Sr) 1 _ 2 .K 2; Fe2As2. 



Ellipsometry has an advantage of measuring the complex dielectric function e(uS) = E\ (oj) +ie2 (u>) = 1 +47tzct(o;) /u) directly 
without a need for reference measurements and Kramers -Kronig transformation (KKT). The KKT allows for a consistency 
check of Ex(oS) and <ji(cj) and implies that Aei at any energy depends on Atri(a;) in the whole spectral range including the 
extrapolation region: 



where £o is the dielectric permittivity, and <Ji(uj) is the real part of optical conductivity. This consistency analysis is rather 
insensitive to the exact shape of the extrapolation chosen but it does fix the total spectral weight 



where luq is the experimental low-energy cutoff frequency. This procedure is illustrated in Fig. 3 c of the main text. Taking the 
experimental difference spectrum of A<7i(oj) with ASW(cl>o) = (blue circle) and carrying out the KKT results in the real part 
of the dielectric function deviating from the measured data (blue line in the inset). Only by increasing the spectral weight below 
ujq = 12 meV to 0.015 eV 2 (red circle) does one achieve complete agreement with the experiment (red line in the inset). The 
exact shape of the extrapolated Aeri (uj) plays a minor role. The maximum uncertainty introduced by the unknown shape can be 
calculated as the difference of two extreme configurations: all ASW(wo) at cj = and u = ljq: 



Kramers-Kronig consistency check 






5Ae£\u) = 8 
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when to ^> ujq. On the other hand, the accuracy to which the spectral weight is determined at the same energy is given by 



*Ae[ a) (w) 



dx 



o+ 



< 



/(ASW(wo)) 



The relative effect of the shape change over the magnitude change of the spectral weight in the extrapolation region is then 
\5Ae < 'i\lo)/SAe < ^\ui)\ — > (luq/lo) 2 , for lo 3> Wo- m tne present case taking ujq = 12 mcV and io = 250 meV (as in the 
inset of Fig. 3 c) one gets a shape uncertainty fraction of 0.2%. Thus the effect is negligible already at rather low frequencies. 
The same analysis applies for the high-energy extrapolation above 6.5 eV. However, complete agreement between Ae\(uj) and 
Aeiiys) within the accuracy of the experiment was found up to 6.5 eV, therefore no experimentally discernable missing spectral 
weight is contained at higher energies. 
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FIG. SI: a Difference spectrum of the real part of optical conductivity (top panel) and dielectric function (bottom panel) between 40 and 10 K 
(red and black open circles, respectively). Lorentzian fit to both spectra (black solid lines). Small background shift (dashed line in the bottom 
panel) obtained from the fit and further confirmed by temperature modulation measurements c-e. b Temperature scan of the real part of the 
optical conductivity between 60 and 15 K at 2.5 eV (blue open circles, right scale) and 24 meV (red filled circles, left scale), c-e Temperature 
modulation of the real part of the optical conductivity (c) and dielectric function (d,e) between 20 K and 40 K with a period of 1800 sec 
averaged over 24 waveforms. 



Temperature-modulation measurements 



The small background shift shown as a dashed line in the bottom panel of Fig. Sit was identified from a simultaneous fit 
of A<7i(u) and Ae\{uj) (upper and bottom panels in Fig. SI i). The superconductivity-induced nature of the suppression of 



the absorption band as well as of the background shift was confirmed by a dynamic temperature scan between 60 and 15 K 



(red open circles in Fig. S1 1). Its temperature dependence clearly follows that of the far-infrared optical conductivity due to 
the opening of the superconducting gap (blue filled circles in Fig. SI :>). To estimate the background shift more accurately, 
temperature modulation measurements of o\ at the resonance photon energy 2.47 eV and e\ at off-resonance photon energies 



of 2.12 eV and 2.82 eV were carried out. In Fig. SI >e the sample temperature was changed between 20 to 40 K with a period 
of 1800 s and later averaged over 24 periods to reduce noise to Aei = 10~ 4 . This confirms a minute background increase of 
A £l (2.47eV) = (8 ± 4) 10~ 4 . 
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LDA calculations: effect of Ef ± A cuts of the density of states on optical conductivity. 



The band-structure calculations were performed using a linear-muffin-tin orbital method in the atomic sphere approxima- 
tion [S3 1 within the LDA starting form the known crystal structure of Ba(Sr) 1 _ 2 ,K K Fe2As2 [Sffl 0. The calculations for 
Ba0.68K0.32Fe2A.S2 predict a contribution to the optical conductivity in the visible that results from transitions with final states 
crossing the Fermi surface, which can experience suppression due to a redistribution of occupied states below the Fermi level 
across the superconducting transition. Apart from the interband transitions shown as two arrows in Fig. 2 c, a contribution from 
transitions with initial states crossing the Fermi surface and into higher-lying hybridized As-p and Fe-d states is significant in 
this energy range. To substantiate that a redistribution of occupied states below the Fermi level within one superconducting 
gap can account for the observed suppression, optical conductivity in the LDA framework was calculated for the two cases of 
hole and electron transitions in the vicinity of T- and M-points of the Brillouin zone, i.e. for the hole pockets of the Fermi 
surface, shown in Fig. S2 a and S2 b, respectively. Shown in the Fig. S2 c is difference optical conductivity spectra between 




FIG. S2: a Difference real part of the optical conductivity between ungapped and gapped regimes for transitions to (red solid line) and from 
(blue dashed line) the Fermi level, b Schematic representation of the optical transitions at T point of the Brillouin zone with a contribution 
to the optical conductivity shown in a. c Population (red) or depletion (blue) of all electronic states within one superconducting gap value 
A = 10 meV. Arrows depict the same optical transitions as in b. 

gapped and ungapped regime. In the gapped case the transitions in the energy window AE = A = 10 meV below (blue 
line) and above (red line) the Fermi level are forbidden, to simulate complete depletion (population) of occupied (empty) states. 
Hole (electron) contribution to the optical conductivity displays suppression only when states above (below) the Fermi level are 
fully populated (depopulated) to within one superconducting gap. The size of the effect in both cases is an order of magnitude 
larger than the experimentally detected suppression. Therefore a fractional population redistribution can indeed account for the 
superconductivity- and SDW-induced anomalies in Bao.6sKo.32Fe2As2 and SrFe2As2. 



Interband optical conductivity in single- and multiband BCS theory 



In the framework of the BCS theory the charge carriers and elementary excitations in the superconducting state differ sig- 
nificantly from those of in the normal state. This is manifested in the modified with respect to free charge carriers dispersion 
of the excitations of the superconducting condensate - Bogolubov quasiparticles. To obtain this dispersion one can introduce 
quasiparticle operators diagonalizing the original BCS Hamiltonian [S|6]|71 

H = J2 £ k™k CT + v ^ c U c -n c -n c n 

ka kl 

as follows: 

c kt = M k7ko + u k7 k i> 

C -14 = - u k7ko+ w k7kr (3) 
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The complex functions of the k-vector u k and w k determine the probability of the pair state comprised of electrons with momenta 
k and — k being empty or occupied, respectively: 




(4) 



where E, 



A 2 is the Bogolubov quasiparticle's dispersion and = e k — jj, is the normal-state electron dispersion 



measured with respect to the chemical potential fi. This probability distributions are inherently smeared around the Fermi level 
in the ground state. This smearing of the quasiparticle occupation probabilities at K resembles closely that of normal-state 
particles at T = T c 

One of the important and quite intuitive consequences of the Bogolubov- Valatin transformation ^ is that the operators 7 k(T 
and c kcr are connected via a unique, one-to-one relation. This immediately implies conservation of the total number of states in 
a given energy range, i.e. 



N sc (E)dE = JV NS (Ode, 



(5) 



where N$c(E) is the quasiparticle density of states (DOS) in the superconducting state, -/Vns(0 i s tne normal-state electron 
DOS. This relation requires that the states within one A below the Fermi level be expelled to energies lower than E-p — A, while 
those within one A above the Fermi level to energies higher than Ep 
conserves the population above and below the Fermi level so that TVgg 



A, as illustrated schematically in Fig. S3 1. This process 



= N^g (blue area equal to the gray area in Fig. |S3p). 



The exact analytical expression for the DOS of Bogolubov quasiparticles follows from equation |5]l bearing in mind that from 
the definition of the quasiparticle dispersion £ k = W £" k — A 2 : 



N S c(E) = JV NS (0 



dE 



(6) 



This expression is plotted in Fig.|S3fr (black solid line) for the case of free-electron normal-state dispersion (gray dashed line). 
The square-root singularity in the quasiparticle DOS stems from the flattening of the normal-state dispersion in the vicinity 
of the Fermi surface, as shown in Fig. S3 3 (the inner hole dispersion corresponds to the DOS plotted in Fig. S3 i). Occupied 



quasiparticle states are depicted as filled black circles, while the quasiparticle vacancies are shown as empty black circles. The 
fading of filled and empty black circles represents the occupation probabilities (|u k | 2 , l^iJ 2 ) in equation J4b. As it has already 
been mentioned above there exists a finite smearing of these probabilities even at K. It leads to finite occupation of those 
regions of the Brillouin zone unoccupied in the normal state, the so-called backfolding of the quasiparticle dispersion, clearly 
visible in Fig. S3 ~>. As optical conductivity only probes averaged over wavevectors regions of the k-space it cannot resolve 



the result of the backfolding as opposed to angle-resolved photoemission spectroscopy (ARPES), where this effect has been 
reliably established [S9-lT|. On the other hand the effect of smearing itself is incorporated in the density of states within 
the BCS formalism and is, therefore, included into our considerations. It is a single-band effect of superconductivity on the 
band structure and does not lead to population redistribution, i.e. the total number of unoccupied states below the transition 



is conserved iVgg 



= iVJJg (blue area is equal to the gray area in Fig. S3 i). This can only lead to a small corrugation of an 



optical absorption band on the scale of one superconducting-gap energy superimposed on the overall broad feature without any 
modification of its spectral weight rS fTZl . The experimentally observed suppression of an absorption band on the scale of its 
full width rather than on the scale of superconductivity-induced modification of the dispersion necessarily requires population 
imbalance ATgg < A^Jg (red area is unequal to the gray area in Fig. S3 1). It is unlikely that such an effect can be a consequence 



of a dynamic population balance between two or more bands at finite temperatures because the temperature dynamics of the 
suppression mimics that of the optical conductivity in the FIR region due to the opening of the superconducting gap, as shown 
inFig.jSTp. 

To account for the multiband character of iron pnictides one may consider the multiband BCS theory [S8 |. It is a straight- 
forward generalization of the single-band BCS theory with the only complication that each separate band has its own gap, 
quasiparticle dispersion, and a pair (u k , u k ). However, in the framework of this multiband theory the quasiparticle operators do 
not involve normal-state particle operators and mix only the creation and annihilation operators from the same band (preferring 
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FIG. S3: (interactive: press buttons to choose) a Density of states in the normal (gray line), conventional superconducting state (black line), and 
an unconventional state with a depletion of unoccupied states (red line). Filled areas of respective colors represent total number of unoccupied 
states, b Schematic representation of the band structure of BKFA in the normal (dashed parabolas) and the superconducting state (filled and 



empty circles for occupied with probability 



and unoccupied with probability 



states with Bogolubov dispersion, respectively). The 



pair (u k 



is unique for each separate band in the conventional multiband BCS approach [S6 8|. The orbitals of dominant contribution to 



each particular band are specified, c (schematic) Difference spectra of the real part of optical conductivity between 40 and 10 K. 



our notation to that of ref. 



--kj. 



'kf 



(c)* , (c) f 
U k e k0 ^k e kl ' 



(c)* 
~ v k e k0 



(c) t 



kl' 



f 4- v (d) f f 



kl- 



,.(<*)* 



f + u {d) f f 
./k0 + "k Jl 



kl' 



(7) 



where (c k ^, d k ^) are normal-state particle operators and (e k ^, / k ^) are the multiband counterparts of the operators 7 k ^ in the 



single-band BCS theory. The coefficients (it 



(c,d) (c,d)> 



certainly depend on the properties of both bands as well as on the inter- 



band coupling but the relations |7]i are still unique one-to-one relations, which immediately implies that, however complicated 
the quasiparticle dispersions may be, for each separate band relation 



(8) 



holds, with (i) running through all bands. As a consequence, just like in the single-band case, the total occupied/unoccupied 
population is conserved across the superconducting transition and thus only changes of interband optical conductivity on the 
scale of 2A are expected. The population imbalance < (red area in Fig. S3 1 smaller than the gray area) required to 



accommodate the experimentally observed suppression of an 2.5 eV absorption band on its full width of about 1 eV can come 
from redistribution of occupation of the different bands below T c . It requires a lowering of the material's chemical potential 
in the superconducting state and, therefore, an additional contribution to the condensation energy. However, the standard BCS 
theory and its generalization to the multiband case do not self-consistently take into account this effect: though predicting a 
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lowering of a chemical potential as a consequence of a non-zero gain in the free energy of the system (condensation energy), 
they premise on equations with an essentially constant chemical potential. Consistent treatment of a variable chemical potential 
might render the Bogolubov-Valatin transformation|7]inappropriate in the multiband case and violate the population conservation 
of the occupied and unoccupied states within each band (as shown in Fig. S3 i) - a fundamental consequence of the standard 
BCS theory. The resulting correction, small as it may be for conventional superconductors, in the presence of large Fe-As bond 
polarizability can lead to a large effect and potentially enhance superconductivity in iron pnictides. 



Spectral weight and kinetic energy gain 

Increase of the condensation energy due to the lowering of the chemical potential of the system in the superconducting state 
explained in the previous chapter can be related to the SW of the suppressed optical band at 2.5 eV. Based on the current 
experimental evidence, it is impossible to confidently specify whether the liberated upon the superconducting transition SW 
contributes to that of the superconducting condensate. Due to the reliably determined background shift of the real part of Aei (cj) 



in the visible spectral range (see Fig. S 1 1, bottom panel) it is certain, however, that this S W is redistributed at energies below 
2 eV. One can easily estimate the order of magnitude of the free-energy gain resulting from this SW transfer assuming that it 
contributes to the intraband itenerant response r ^T3l[T4ll : ASW(Q.) = (ne 2 a 2 /2h 2 V u ) (-E K ), where SW(ft) = J 0+ ai(ui)du) 
is the in-plane spectral weight, a is the in-plane lattice constant, V u - unit cell volume and (— £k) — 0-/N) J^k a n k J^^l^^a 
is a measure of the system's kinetic energy per unit cell. In the case of the tight-binding nearest-neighbor approximation 
it is exactly equal to the kinetic energy of the charge carriers per unit cell and thus contributes directly to the condensation 
energy. This approximation might not hold for iron-pnictide superconductors but the above estimate shows that though small 
this additional spectral weight in the superconducting Bao.68Ko.32Fe 2 As 2 is enough to account for the condensation energy in 
this compound. 



* Correspondence should be addressed to A. C. and A. V. B.: a.charnukha@fkf.mpg.de, a.boris@fkf.mpg.de 
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